Wind morphology around cool evolved stars in binaries: the case of slowly accelerating oxygen-rich outflows

The late stellar evolutionary phases of low and intermediate-mass stars are strongly constrained by their mass-loss rates. The wind surrounding cool evolved stars frequently shows non-spherical features, thought to be due to an unseen companion orbiting the donor star. We study the morphology of the circumbinary envelope, in particular around oxygen-rich asymptotic giant branch (AGB) stars. We run a grid of 70 3D hydrodynamics simulations of a progressively accelerating wind propagating in the Roche potential formed by a mass-loosing evolved star in orbit with a main sequence companion. We resolve the flow structure both in the immediate vicinity of the secondary, where bow shocks, outflows and wind-captured disks form, and up to 40 orbital separations, where spiral arms, arcs and equatorial density enhancements develop. When the companion is deeply engulfed in the wind, the lower terminal wind speeds and more progressive wind acceleration around oxygen-rich AGB stars make them more prone than carbon-rich AGB stars to display more disturbed outflows, a disk-like structure around the companion and a wind concentrated in the orbital plane. In these configurations, a large fraction of the wind is captured by the companion which leads to a significant shrinking of the orbit over the mass-loss timescale, if the donor star is at least a few times more massive than its companion. Provided the companion has a mass of at least a tenth of the mass of the donor star, it can compress the wind in the orbital plane up to large distances. Our grid of models covers a wide scope of configurations function of the dust chemical content, the terminal wind speed relative to the orbital speed, the extension of the dust condensation region around the cool evolved star and the mass ratio. It provides a frame of reference to interpret high-resolution maps of the outflows surrounding cool evolved stars

SPH modelling of companion-perturbed AGB outflows including a new morphology classification scheme

CONTEXT: Asymptotic giant branch (AGB) stars are known to lose a significant amount of mass by a stellar wind, which controls the remainder of their stellar lifetime. High angular-resolution observations show that the winds of these cool stars typically exhibit mid- to small-scale density perturbations such as spirals and arcs, believed to be caused by the gravitational interaction with a (sub-)stellar companion. AIMS: We aim to explore the effects of the wind-companion interaction on the 3D density and velocity distribution of the wind, as a function of three key parameters: wind velocity, binary separation and companion mass. For the first time, we compare the impact on the outflow of a planetary companion to that of a stellar companion. We intend to devise a morphology classification scheme based on a singular parameter. METHODS: We ran a small grid of high-resolution polytropic models with the smoothed particle hydrodynamics (SPH) numerical code PHANTOM to examine the 3D density structure of the AGB outflow in the orbital and meridional plane and around the poles. By constructing a basic toy model of the gravitational acceleration due to the companion, we analysed the terminal velocity reached by the outflow in the simulations. RESULTS: We find that models with a stellar companion, large binary separation and high wind speed obtain a wind morphology in the orbital plane consisting of a single spiral structure, of which the two edges diverge due to a velocity dispersion caused by the gravitational slingshot mechanism. In the meridional plane the spiral manifests itself as concentric arcs, reaching all latitudes. When lowering the wind velocity and/or the binary separation, the morphology becomes more complex: in the orbital plane a double spiral arises, which is irregular for the closest systems, and the wind material gets focussed towards the orbital plane, with the formation of an equatorial density enhancement (EDE) as a consequence. Lowering the companion mass from a stellar to a planetary mass, reduces the formation of density perturbations significantly. CONCLUSIONS: With this grid of models we cover the prominent morphology changes in a companion-perturbed AGB outflow: slow winds with a close, massive binary companion show a more complex morphology. Additionally, we prove that massive planets are able to significantly impact the density structure of an AGB wind. We find that the interaction with a companion affects the terminal velocity of the wind, which can be explained by the gravitational slingshot mechanism. We distinguish between two types of wind focussing to the orbital plane resulting from distinct mechanisms: global flattening of the outflow as a result of the AGB star’s orbital motion and the formation of an EDE as a consequence of the companion’s gravitational pull. We investigate different morphology classification schemes and uncover that the ratio of the gravitational potential energy density of the companion to the kinetic energy density of the AGB outflow yields a robust classification parameter for the models presented in this paper

SPH modelling of wind-companion interactions in eccentric AGB binary systems

The late evolutionary stages of low- and intermediate-mass stars are characterised by mass loss through a dust-driven stellar wind. Recent observations reveal complex structures within these winds, that are believed to be formed primarily via interaction with a companion. How these complexities arise, and which structures are formed in which type of systems, is still poorly understood. Particularly, there is a lack of studies investigating the structure formation in eccentric systems. We aim to improve our understanding of the wind morphology of eccentric AGB binary systems by investigating the mechanism responsible for the different small-scale structures and global morphologies that arise in a polytropic wind with different velocities. Using the smoothed particle hydrodynamics (SPH) code Phantom, we generate nine different high-resolution, 3D simulations of an AGB star with a solar-mass companion with various wind velocity and eccentricity combinations. The models assume a polytropic gas, with no additional cooling. We conclude that for models with a high wind velocity, the short interaction with the companion results in a regular spiral morphology, that is flattened. In the case of a lower wind velocity, the stronger interaction results in the formation of a high-energy region and bow-shock structure that can shape the wind into an irregular morphology if instabilities arise. High-eccentricity models show a complex, phase-dependent interaction leading to wind structures that are irregular in three dimensions. However, the significant interaction with the companion compresses matter into an equatorial density enhancement, irrespective of eccentricity.Comment: 23 pages, 22 figure

MAGRITTE, a modern software library for 3D radiative transfer: I. Non-LTE atomic and molecular line modelling

Radiative transfer is a key component in almost all astrophysical and cosmological simulations. We present MAGRITTE: a modern open-source software library for 3D radiative transfer. It uses a deterministic ray-tracer and formal solver, i.e. it computes the radiation field by tracing rays through the model and solving the radiative transfer equation in its second-order form along a fixed set of rays originating from each point. MAGRITTE can handle structured and unstructured input meshes, as well as smoothed-particle hydrodynamics (SPH) particle data. In this first paper, we describe the numerical implementation, semi-analytic tests and cross-code benchmarks for the non-LTE line radiative transfer module of MAGRITTE. This module uses the radiative transfer solver to self-consistently determine the populations of the quantized energy levels of atoms and molecules using an accelerated Lambda iteration (ALI) scheme. We compare MAGRITTE with the established radiative transfer solvers RATRAN (1D) and LIME (3D) on the van Zadelhoff benchmark and present a first application to a simple Keplerian disc model. Comparing with LIME, we conclude that MAGRITTE produces more accurate and more precise results, especially at high optical depth, and that it is faster

Relativistic Beaming and the Intrinsic Properties of Extragalactic Radio Jets

Relations between the observed quantities for a beamed radio jet, apparent transverse speed and apparent luminosity (beta_app,L), and the intrinsic quantities, Lorentz factor and intrinsic luminosity (gamma,L_o), are investigated. The inversion from measured to intrinsic values is not unique, but approximate limits to gamma and L_o can be found using probability arguments. Roughly half the sources in a flux density--limited, beamed sample have a value of gamma close to the measured beta_app. The methods are applied to observations of 119 AGN jets made with the VLBA at 15 GHz during 1994-2002. The results strongly support the common relativistic beam model for an extragalactic radio jet. The (beta_app,L) data are closely bounded by a theoretical envelope, an aspect curve for gamma=32, L_o= 10^25 W/Hz. This gives limits to the maximum values of gamma and L_o in the sample: gamma_max about 32, and L_o,max ~ 10^26 W/Hz. No sources with both high beta_app and low L are observed. This is not the result of selection effects due to the observing limits, which are flux density S>0.5 Jy, and angular velocity mu<4 mas/yr. Many of the fastest quasars have a pattern Lorentz factor gamma_p close to that of the beam, gamma_b, but some of the slow quasars must have gamma_p<<gamma_b. Three of the 10 galaxies in the sample have a superluminal feature, with speeds up to beta_app about 6. The others are at most mildly relativistic. The galaxies are not off-axis versions of the powerful quasars, but Cygnus A might be an exception.Comment: 12 pages, 9 figures, 1 table, accepted for publication in the Astrophysical Journa

Prazosin during threat discrimination boosts memory of the safe stimulus

The α-1 adrenoreceptor antagonist prazosin has shown promise in the treatment of post-traumatic stress disorder (PTSD) symptoms, but its mechanisms are not well understood. Here we administered prazosin or placebo prior to threat conditioning (day 1) and tested subsequent extinction (day 2) and reextinction (day 3) in healthy human participants. Prazosin did not affect threat conditioning but augmented stimulus discrimination during extinction and reextinction, via lower responding to the safe stimulus. These results suggest that prazosin during threat acquisition may have influenced encoding or consolidation of safety processing in particular, subsequently leading to enhanced discrimination between the safe and threatening stimuli

MAGRITTE, a modern software library for 3D radiative transfer - II. Adaptive ray-tracing, mesh construction, and reduction

Radiative transfer is a notoriously difficult and computationally demanding problem. Yet, it is an indispensable ingredient in nearly all astrophysical and cosmological simulations. Choosing an appropriate discretization scheme is a crucial part of the simulation, since it not only determines the direct memory cost of the model but also largely determines the computational cost and the achievable accuracy. In this paper, we show how an appropriate choice of directional discretization scheme as well as spatial model mesh can help alleviate the computational cost, while largely retaining the accuracy. First, we discuss the adaptive ray-tracing scheme implemented in our 3D radiative transfer library MAGRITTE, that adapts the rays to the spatial mesh and uses a hierarchical directional discretization based on HEALPIX. Second, we demonstrate how the free and open-source software library GMSH can be used to generate high-quality meshes that can be easily tailored for MAGRITTE. In particular, we show how the local element size distribution of the mesh can be used to optimize the sampling of both analytically and numerically defined models. Furthermore, we show that when using the output of hydrodynamics simulations as input for a radiative transfer simulation, the number of elements in the input model can often be reduced by an order of magnitude, without significant loss of accuracy in the radiation field. We demonstrate this for two models based on a hierarchical octree mesh resulting from adaptive mesh refinement, as well as two models based on smoothed particle hydrodynamics data